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Abstract 

We present the exact numerical solutions for the Schwinger-Dyson equations at finite temper¬ 
ature with general gauge in Abelian gauge theory. We then study the chiral phase transition on 
temperature from the obtained solutions. We find that, within the quenched Schwinger-Dyson 
equations, there exists substantial gauge dependence on the solutions and the critical temperature. 
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I. INTRODUCTION 


The chiral phase transition at finite temperature is an important and interesting phe¬ 
nomenon in gauge theories. The chiral symmetry is broken for low temperature in the 
strongly coupling system, and expected to be restored for high temperature. Then, it is 
interesting to investigate this critical behavior; in particular, the critical temperature on the 
chiral phase transition. 


The Schwinger-Dyson equation (SDE) |l| is frequently used method for the studies on the 
above mentioned chiral phase transition. The SDE are the integral equations for the green 


functions, and the solutions lead t 


generated mass (for reviews, see, 
on the SDE at zero temperature 


le field strength renormalization and the dynamically 


H 


. A lot of works have been devoted to the analyses 


and finite temperature 


14 


2l| . The investigations 


at finite temperature are rather limited comparing to the ones at zero temperature, since it 
becomes difficult to solve the equations at finite temperature due to the separation of the 
time and space components 22(]. Then the preceding works have mostly been done through 
employing some approximations on the equations in a specific gauge. However, it is possible 


to evaluate the equations without approximations 


14| . In this paper, we shall solve the 


SDE at finite temperature without approximations in general gauge. Thereafter we study 
the gauge dependence on the chiral phase transition with respect to temperature. 

The plan of the paper is as follows. Section Ullpresents the equations at finite temperature 
and the prescription of solving the equations. The numerical solutions and the critical 
behavior are shown in Secs. im and IIVI The concluding remarks are given in Sec. |Vl 


II. SCHWINGER-DYSON EQUATION 

We will here present the SDE at finite temperature and the iteration method which is 
the numerical procedure employed in this letter. 


A. The equations at finite temperature 

The equation for the fermion self-energy S(p) is given by 

S(p)=*e"/(1) 
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with the coupling strength, e, the propagators for the photon and fermion, D^^{p — q) and 
S{q), and the vertex function r'^(p, q). For and F*^, we use the following quenched forms 


=-^ + (1 - 0^ 

T'^{p,q) = Y, 


( 2 ) 

(3) 


where ^ is the gauge parameter and = P^i — Qfi- 
The fermion propagator is defined by 

S{po,p) = 


C{po,phoPo + A(po,p) 7 iP* - B{po,p) 


(4) 


with p = |p|. Note that the time and space components are to be treated separately for the 
extension to hnite temperature. In the imaginary time formalism, the continuous qn integral 
is replaced by the discretized summation due to the applied boundary condition |22l |. 


1 /■°° 

— / dqoF{qo) T 
2vr* 


OO 

E 

m=—oo 


F (ibJm) 


(5) 


where T is the temperature of the system. The frequency, Um, is taken as 


ujm = (2m + l)7rT 


( 6 ) 


for the fermionic held. 

After some algebras one obtains the following equations 


C{0Jn,p) = 1 + -T / dg [XccC{uJm, q) + XcAA{oJra, 

^ J\ 

m=—oo ^ 

^ OO /-A 

A{Un,p) = 1 + -T ^ / dg [XAcC{Um,q) +XAAA{uJm, 

TT — \ 

m=—oo '' 

B{un,p) = mo +-T / dq [XBB{um,q)]A{u^,q), 

^ Jx 


m=—oo 


with a = ! (dvr). 


A(a;m,g) = 


C 2 (a;m, g)a;^ + A^{um, q)q^ + B'^{um, 


g)] A(a;m,g), 

(7) 

g)] A(a;,n,g), 

(8) 


(9) 

q)' 

(10) 
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and 


Xcc = + 

OJri. OJrt. 


XcA = e-—/i + e-— 


OJr, 


UJrt 


2 




q +p" 


h, 


2 (- ^ m T A 

AC — - 4 -——h - 4 - 

Xaa = 




p2 ^ p2 

2 g 2 1 


uj'J + q^-p^ 


^3 ‘^'m + ^+( 5 '^ + P^) h - ■^^- 


p 2 2 p 2 

Is = -53 +/i, 

2 p uj'J + {q + pf 


/ 4 / 2 2 n 2 

- (9 - P ) 


/ 2 , 


12 = 


q 


2p [u'J + (g - pf U}'J + {q + pf_ ■ 


( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 


Here = a;„ — cUm, 4± = 1 ± 4 and 4^ = 3 ± 4, and we introdnced the infrared and 
ultraviolet cntoffs, A and A. These eqnations are the set of the Schwinger-Dyson eqnations 
at hnite temperatnre which are the straightforward extension from the eqnations at zero 
temperatnre with three dimensional momentnm cntoff 23 1 . 

Once the solntions are obtained, we can evalnate the dynamical mass and the chiral 
condensate throngh the relations, 


M{Un,p) = B{Un,p)/A{Un,p), 

/If, 

m — — rv^ '' V / 


tr 


B{uim, q) 


C^Cm, q)iqm + q)q^ + q) 


(18) 

(19) 


It shonld be noted that the chiral condensate becomes zero when the dynamical mass van¬ 
ishes. Therefore one can see whether the chiral symmetry is broken, 0 7 ^ 0, from the valne of 
M{un,p)- Strictly, the phase transition should be evaluated at the minimum of the effective 
potential which is derived by integrating out the gap equation with respect to the order 


parameter. However, it is known to be ac 


equate to see the solution of the gap equation in 


the system at zero chemical potential 17|, then we search the critical temperature on the 


chiral phase transition by checking the value of the dynamical mass in this paper. 
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B. The iteration method 

In this subsection, we discuss the numerical procedure for the integral equations. To 
briefly present the basic idea of the iteration method, we show the equations with the 
trapezoidal rule. (In the actual calculations, we used the Gauss-Legendre integration for the 
g—integration.) 

In solving the equations, we discretize the integral to discrete summation as 

Nt-l Np-1 

_ 1 r [t'C'C' Ml-^) \jCA A^-l) I'on'i 

^nj ^ / j ^ / j ^njmk^mk ' -^njmk^mk ^mk i 

rn=—Nt k=0 

Nt-l Np-l 

^ / J ^ , ^njmk^mk ‘ -^njmk^mk ^mk ’ G-*-/ 

m=—Nt k=0 
Nt-l Np-1 

B™ = mo + ^11 (22) 

m=—Nt k=0 

with (5 = (A — X)/Np. Note that we introduce the cutoff on the frequency summation since 
it is technically impossible to take the inhnite number of summation. However, as will be 
conhrmed by the numerical calculations in the next section, the contributions from the large 
number of n and m are negligible, then the artifact of setting Nt does not affect the solutions 
if we choose large enough number for Nt- The matrix of the form Fnjmk indicates 

Fnjmk F{}^ni Pi ^mi 0^ 1 P A T 6j, Q A T 6k, (^3) 

and the superscript / in F^’'^ means the number of iteration. Therefore, the / th solutions, 
are derived from the equations with the previous solutions, F^’'~^\ By setting the initial 
conditions on and then carrying out the enough times of iterations, one can 

obtain the solutions. 

III. NUMERICAL SOLUTION 

Having presented the equations and the procedure of solving them, we are now ready to 
perform the numerical analyses. In this section, we are going to show the numerical solutions 
for C, A and B with respect to Un, p, and temperature. 
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A. Solutions for C{un,p), A{u}n,p), B{ujn,p) 


Figure [T] displays the solutions of C{un,p), A{un,p), B{un,p) for a = 3.5, A = O.OIA, 
mo = 0 with ^ = 3, 1 and 0 at temperatures T = 0.05A and 0.15A. In carrying out the 
calculations, we took the number of frequency summation Nt = 50 and performed the 100 
times of iterations. We have numerically conhrmed that these numbers are good enough in 
our present analyses, namely the solutions well converge with W = 50 and W = 100. Then, 
we will persistently use this setting in what follows. 

From Fig. [H one sees that C, A and B decrease with increasing temperature, and these 
shapes become smoother when temperature increases. This result on temperature can be 
understood easily, since the non-perturbative effect is expected to be suppressed at high 
temperature. Therefore the symmetry tends to be restored at high temperature, which can 
be conhrmed by the obtained solutions. On the gauge dependence, the values of C and 
A are large for large while B becomes smaller with respect to It is known from the 
T = 0 analyses with the four dimensional momentum cutoff regularization that C(= A) and 
A become 1 in the Landau gauge = 0. Then the deviation from 1 is expected to become 
larger when ^ increases, which is conhrmed by the obtained results. On the other hand, the 
value of B is small for large This comes from the fact that the denominator of the fermion 
propagator becomes larger when C and A raise up, then the integral in the equation for B 
becomes small, leading smaller value for B consequently. 


B. Temperature dependence on the solutions 

For the purpose to see the temperature dependence, it may be nice to display the hgures 
with the temperature axis. Figures [2] and [3] show the temperature dependence on C{un, A), 
A(a;„,A), B{un,X), and C{uo,p), A{uo,p), B{uo,p). We observe that C and A slightly 
increases according to T when B 0, then decreases after B vanishes. While the result of 
B is easier to see; it stays almost constant for low T, then drops at certain T. This is the 
clear signal of the phase transition, where the mass factor B becomes zero at some critical 
temperature T^. The issue is interesting; we will discuss on this in more detail in the next 
section. 
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C(a)„,p) C((fl„,p) C(o)„,p) C((o„,p) C((o„,p) C((fl„,p) 


^=3, T=0.05A 




















FIG. 1. C{ujn,p), A{ujn,p), B{(jJn,p) with a = 3.5, for T = 0.05A and T = 0.15A. 
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C((Oo,p)/A C((Oo,p)/A C((Oo,p)/A C(co„, X) C(co„, X) C(co„, X) 



FIG. 3. C{uJo,p), A{u;o,p), B{uo,p) for a = 3.5. 
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IV. CHIRAL PHASE TRANSITION 


We have shown rough shapes of the solutions by aligning the three dimensional hgures 
in the previous section. We will now discuss the gauge dependence on the chiral phase 
transition in more detail focusing on the solutions at the lowest frequency and momentum, 
i.e., the values at {ujn,p) = A). 


A. Gauge dependence on the solutions 

Figure m shows how the solutions C{u)o, A), A{ojq, A), B{ojo, A) and M{u)o, A) changes with 
respect to T. As we previously saw that C increases with T when B is nonzero, then 






FIG. 4. C, A, B and M as the functions of T with a = 3.5 for ^ = 0, 1 and 3 . 

decreases after B vanishes. A shows the similar curves for = 0 and 1, while it almost 
monotonously decreases for ^ = 3. B and M become small for high T as already observed 
in Figs. [2]and[3l 

From the hgure, one can clearly conhrm the quantitative gauge dependence on the solu¬ 
tions. C and A are around 3 for low T with = 3, and within the range of around 1 — 2 
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for ^ = 0 and 1. M at T = 0 is around 0.7 for ^ = 0 and it becomes considerably smaller 
for = 3. Thus, there exists the indispensable gauge dependence on the solutions for the 
quenched SDE. 


B. Critical temperature 


Let us study the critical temperature, T^, of the chiral phase transition here. We have 
discussed that the dynamically generated mass can be regarded as the order parameter of 
the transition in Sec. Ill Al Then we shall carefully search the critical temperature of the 
chiral phase transition by evaluating the value of M. 






FIG. 5. M{ujq, A) as the functions of T for various ^ and a. 

Figure [5] displays the temperature dependence on the dynamical mass, M{uo,X), for 
various ^ and a. The critical temperature becomes larger when a increases. This can be 
easily understood since the system has stronger correlation for larger coupling. While is 
smaller for larger This comes from the effect of the held renormalization as we have seen 
from the preceding section. The explicit numbers of the critical temperature for various ^ 
and a are aligned in Tab. [B 
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TABLE I. Critical temperature (Tc/A). 



a = 2 

a = 3 

0 = 4 

0 = 5 

o 

II 

Uy’ 

0.089 

0.187 

0.267 

0.341 

e = i 

— 

0.126 

0.203 

0.265 

II 

— 

0.07 

0.162 

0.223 

II 

CO 

— 

— 

0.14 

0.197 


Finally, we show as the function of a in Fig. [ 6 l We again conhrm that there exists 



FIG. 6. Critical temperature with respect to a. 


substantial gauge dependence on the critical behavior. Note that we can not extend the 
curves for very low temperature within this formalism, since we need to take large number 


of frequency summation to get reliable results for low 


24| . which requires enormous 


calculation time. However one easily sees the limiting values on the critical coupling for 
Tc —>■ 0 from the hgure, which are around Oc — 1.2, 2, 2.6 and 3 for ^ = 0, 1, 2 and 3. It 
may be interesting to mention that these values are close to the ones obtained in the SDE 
with the four dimensional cutoff regularization at T = 0, which leads Oc — vr/3, 2 and 3 for 


.^ = 0, 1 and 3 


13| 


V. CONCLUDING REMARKS 

We have evaluated the numerical solutions for the SDE at hnite temperature without 
approximations in Abelian gauge theory in this paper. There we observe the considerable 
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gauge dependence on the solutions. We also analyzed the critical behavior of the chiral 
phase transition within the SDE approach. We conhrm that the critical temperature as well 
drastically depends on the chosen gauge. This indicates that, when one investigates the 
chiral phase transition in the quenched SDE with using some approximations, one should 
be very careful on the set of assumptions. In particular, the field renormalization factors C 
and A nonnegligibly deviate from the unity for nonzero 

Since we started from the SDE with the quenched form, then the gauge dependence 
as observed in this analysis is inevitable. This is clearly unsatisfactory feature of the 
quenched equations. Therefore the further analyses with generalized equations, such as 


the unquenched investigations 


13|, or the extended vertex [l^, ll2| and coupling 


0 . 


are 


obviously needed; this should be the future direction in the analyses based on the Schwinger- 
Dyson equation. 
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